The self-injection threshold in self-guided laser wakefield accelerators 
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A laser pulse traveling through a plasma can excite large amplitude plasma waves that can be used 
to accelerate relativistic electron beams in a very short distance — a technique called laser wakefield 
acceleration. Many wakefield acceleration experiments rely on the process of wavebreaking, or self- 
injection, to inject electrons into the wave, while other injection techniques rely on operation without 
self-injection. We present an experimental study into the parameters, including the pulse energy, 
focal spot quality and pulse power, that determine whether or not a wakefield accelerator will self- 
inject. By taking into account the processes of self-focusing and pulse compression we are able to 
extend a previously described theoretical model, where the minimum bubble size k p rb required for 
trapping is not constant but varies slowly with density and find excellent agreement with this model. 



Laser wakefield acceleration, where an intense laser 
pulse drives a plasma wave with a relativistic phase ve- 
locity, is a promising technique for the development of 
compact, or "table-top", particle accelerators and radia- 
tion sources. Plasma waves driven in moderate density 
plasmas can support electric fields over a thousand times 
stronger than those in conventional accelerators. Laser 
driven plasma waves have demonstrated electron accel- 
eration to ^ 1 GeV in distances ~ 1 cm [T]-[3]. These 
compact particle accelerators have significant potential as 
bright x-ray sources [3H6] offering peak brightness com- 
parable to 3rd generation synchrotron sources in x-ray 
flashes on the order of just 10 fs. 

At the heart of the laser wakefield acceleration con- 
cept is the fact that electron plasma waves with relativis- 
tic phase velocities are driven to very large amplitudes, 
where they become highly non-linear. If the plasma 
wave is driven beyond a threshold amplitude, the wave 
breaks. When the wave is driven far beyond the wave- 
breaking threshold, the wave structure is destroyed and 
large amounts of charge can be accelerated to high energy 
but with a broad energy spread [7]. With appropriately 
shaped laser pulses this normally catastrophic process 
of wavebreaking can be tamed to produce high quality 
beams of electrons. This is because close to the wave- 
breaking threshold the nature of wavebreaking changes 
- some electrons from the background plasma can be- 
come trapped in the wave without destroying the wave 
structure, a process called self-injection. 

The highly non-linear broken wave regime [8] is used in 
many experiments to produce quasi-monoenergetic elec- 
tron beams [9HTT]. In such experiments a threshold 
plasma density is commonly observed, below which no 
electron beams are produced. Due to the inverse scal- 
ing of the electron beam energy with plasma density, the 
highest energy beams achievable with a given laser sys- 
tem are achieved just above the threshold, and it is well 



known that many of the beam parameters including the 
spectrum and stability are also optimised just above the 
threshold density [I2j [T3] , It is also well known that 
to achieve self-injection at lower densities higher power 
lasers are required - although the exact scaling of the 
threshold with laser power is not well known. A number 
of techniques to improve the electron beam parameters 
including stability and total charge, have recently been 
demonstrated by using alternative injection schemes [TH - 
US] . Crucially these schemes all rely on operating the 
LWFA below the self-injection threshold. A number of 
recent purely theoretical papers have addressed the dy- 
namics of wavebreaking or self-injection [T9H22] . Clearly 
a good understanding of the self-injection threshold is 
important for the development of laser wakefield acceler- 
ators. We report here on a series of experiments which 
identify the key laser and plasma parameters needed to 
predict the density threshold and we develop a model ca- 
pable of predicting the self-injection threshold density for 
a given set of experimental parameters. 

In LWFA experiments the laser pulse self-focuses due 
to the transverse non-linear refractive index gradient of 
the plasma [23] [24] and the spot size decreases towards a 
matched spot size. This matched spot size occurs when 
the ponderomotive force of the laser balances the space 
charge force of the plasma bubble formed. In situations 
where there is no loss of energy during self-focusing, nor 
any change in the pulse duration, the final matched spot 
size, and hence the final intensity is simply a function 
of aP/P c . P is the laser power; a is the fraction of 
laser energy within the full width at half maximum in- 
tensity of the focal spot — important because energy in 
the wings of the spot are not self-focused by the plasma 
wave and so do not contribute; P c is the laser power 
where relativistic self-focussing dominates over diffrac- 
tion, P c = (87re ra 2 c 5 /e 2 )(n c /n e ) ~ 17n c /n e GW (where 
n e is the background plasma electron density and n c 
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is the critical density for propagation of the laser in 
the plasma). We might therefore expect that the self- 
injection threshold would occur at a fixed value of aP/P c 
[25] . However it is also known that the longitudinal non- 
linear refractive index gradient also has a significant ef- 
fect on the pulse properties [26, 27 and we expect this 
to have an affect on the self-injection threshold. 

The experiment was carried out using the multi-TW 
laser at the Lund Laser Centre. The laser delivered pulse 
energies of up to 0.7 J in pulses as short as 40 fs, cor- 
responding to a peak power of 18 TW. An f/9 off- axis 
parabolic mirror was used to focus the pulse. A de- 
formable mirror was used to optimise the focal spot, pro- 
ducing a spot size of 16 ± 1 fim fwhm. For a gaussian 
focal spot the theoretical maximum fraction of energy 
within the fwhm is a — 1/2, the best focus that we ob- 
tained had a = 0.48. The focal plane was positioned 
onto the front edge of a supersonic helium gas jet with 
an approximately flat-top profile of length 1.8 ± 0.1 mm. 

To investigate the self-injection threshold we studied 
the effect of the plasma density, n e , the total laser energy, 
E, the focal spot quality, a, and the pulse duration, r, on 
the amount of charge in the electron beam. We chose to 
use the total charge in the electron beam as the diagnostic 
of self-injection as it provides a clear unambiguous signal 
of an electron beam. 

The charge was measured using an electron beam pro- 
file monitor, consisting of a lanex screen placed on the 
back surface of a wedge (which was used to collect the 
transmitted laser light). The wedge was 1 cm thick and 
made of glass and therefore prevented electrons below ap- 
proximately 4 MeV reaching the lanex. The lanex screen 
was imaged onto a 12 bit CCD camera. To reduce the 
amount of background light from the interaction, a nar- 
row band interference filter matched to the peak emission 
of the lanex screen was placed in front of the camera. In 
addition the camera was triggered several microseconds 
after the interaction but within the lifetime of the lanex 
fluorescence. The lanex screen was calibrated using the 
absolute efficiency data, absolute response of the CCD 
camera and the details of the imaging system [28] A beam 
profile monitor was used in preference to an electron spec- 
trometer due to the fact that it has a higher sensitivity 
(i.e. the signal produced by a low charge beam dispersed 
inside a spectrometer will drop below the background 
level, whereas the same low charge beam will produce a 
bright image on the profile monitor). Also close to the 
threshold we do not expect the electrons to have partic- 
ularly high energy (i.e. injection could be occurring but 
the electron beam energy could be outside the range of 
the electron spectrometer). . 

The gas jet could produce electron densities up to 
n e = 5 x 10 19 cm -3 . The laser pulse energy was varied 
by altering the energy pumping the final laser amplifier. 
We used the deformable mirror to reduce a by adding 
varying amounts of spherical aberration. Spherical aber- 
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FIG. 1. (colour online) Electron beam profiles for various 
plasma densities for different values of the amount of laser 
energy within the fwhm of the focal spot, a), b) and d) 
kept the total laser energy constant but varied a whereas c) 
reduced the laser energy. Each panel is an average of 5 shots 
and is displayed on a logarithmic colourscale 



ration has the effect of decreasing a without introducing 
asymmetry to the focal spot and without significantly af- 
fecting its size. Degrading the focal spot symmetrically 
was desirable as asymmetric pulses can drive asymmet- 
ric wakes which can have a strong effect on the dynamics 
of self-injection [29]. The pulse duration was altered by 
changing the separation of the gratings in the compres- 
sor. Changing the grating separation introduced both a 
chirp to the pulse spectrum and a skew to the pulse en- 
velope. To take this into account we investigated both 
positive and negative chirps. 

Figure [I] shows the effect of varying the laser pulse en- 
ergy within the focal spot on the self-injection threshold. 
Keeping the total laser energy constant and degrading 
the focal spot (i.e. lowering a) moves the threshold to 
higher plasma densities. We also observe an increase in 
the threshold density when we keep a constant and re- 
duce the laser pulse energy. In fact we find that the two 
effects are equivalent, i.e. that the threshold shifts ac- 
cording to the product aE. This demonstrates that it is 
only the energy within the FWHM of the focal spot that 
contributes to driving the plasma wave. This emphasises 
the importance of laser focal spot quality in LWFA ex- 
periments [30], which are often performed with a ~ 0.3 
[2] [25] . Improving the focal spot could therefore result in 
a significant increase in the electron beam energy achiev- 
able from a given laser system. 

The observed variation of the threshold with aE is 
as expected for one based on aP/P c but this can only 
be confirmed by the behaviour of the threshold when 
we vary the laser pulse duration, keeping aE constant. 
When we do this we see markedly different behaviour. 
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We kept the plasma density constant, at a value just 
above the threshold density for the optimally compressed 
pulse. At this density (n e = 1.6 x 10 19 cm -3 ), with full 
laser energy (aE = 0.32 J) and the fully compressed 
pulse (r = 42 fs) we observed a bright electron beam. 
When we reduced either the plasma density or the pulse 
energy by a small factor (20 - 25%) this beam disap- 
peared, i.e. we dropped below the threshold. Even after 
increasing the pulse duration by a factor of two electrons 
are clearly still injected, as shown in fig [2] This is true 
regardless of the chirp of the laser pulse, however we do 
see an enhancement of the total charge using positively 
chirped (red at the front) pulses as reported previously 
[31]. These pulses have a fast rising edge indicating that 
the precise shape of the pulse may play a role in the to- 
tal charge injected. The direction of chirp of the pulse 
may also affect the rate at which pulse compression oc- 
curs [32 . For both directions of chirp the fact that the 
threshold behaviour is so significantly different to that 
observed when varying aE suggests that pulse compres- 
sion is indeed playing an important role in determining 
whether or not the accelerator reaches wavebreaking. 

In figure [3] we plot the total charge observed on the 
profile monitor screen for the various data sets. Fig- 
ure^) shows the total charge, plotted against the pulse 
power normalised to the critical power for self-focusing, 
for the data sets where we varied the plasma density and 
the energy within the focal spot (either by varying the 
spot quality a or total pulse energy E). The charge rises 
rapidly with increasing aP/P c until eventually reaching 
a plateau at around aP/P c ~ 4. There is an increase in 
the total charge of a factor of ten between aP/P c = 2 
and aP/P c = 4 for both sets of data. The fact that both 
datasets lie on the same curve confirms the fact that it 
is the energy within the focal spot which determines the 
wakefield behaviour. This supports the hypothesis that 
energy in the wings of the focal spot is not coupled into 
the accelerator: energy in the wings of the spot is effec- 
tively wasted. 

Figure [3)3) shows the charge plotted against aP/P c for 
a data set where we kept the plasma density and aE 
constant but varied the pulse duration (by introducing 
either positive or negative chirp) . The markedly different 
behaviour is once again apparent: rather than the rapid 
increase of charge between aP/P c = 2 and aP/P c = 4 
the charge is approximately constant for each data set. 

Figure]^) plots all of the data sets (varying a, E and 
r) against a scaled pulse energy aEn e /n c rather than 
the scaled pulse power. The fact that the pulse dura- 
tion dataset now fits closely with the aE datasets con- 
firms that pulse compression is playing an important role 
in determining whether or not the wakefield accelerator 
reaches self-injection. 

A recent paper that examined the trajectory of elec- 
trons inside the plasma bubble [21] predicts that self- 
trapping will occur when the radius of the plasma bubble 
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FIG. 2. Electron beam profiles for various pulse durations at 
fixed aE and at a plasma density just above the threshold 
density for injection for 40 fs pulses. The pulse duration was 
varied by changing the compressor grating separation which 
introduces a chirp to the pulse a) negative chirp b) positive 
chirp 



(r*b) is larger than a certain value given by: 

k p r b > 2^1n(2 7 |) - 1 (1) 



Where j p ~ \/n c / (3n e ) [33] is the Lorentz factor asso- 
ciated with the phase velocity of the bubble. When this 
condition is met, an electron starting at rest a distance 
from the laser axis and following an elliptical trajectory 
in the bubble fields (thus defining the edge of the bubble) 
will be accelerated by the bubble fields up to 7 P m e c 2 by 
the time it reaches the back of the bubble. A key feature 
of this model is that the normalised bubble size required 
for self-injection k p rb is not constant with density. As 
equation [l] depends only on the plasma density and bub- 
ble size we can determine the minimum pulse properties 
required to reach the threshold by noting that the radius 
of the bubble is related to the pulse energy and duration 
through [34] : 
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Combining equations [T] and [2] yields an expression for the 
minimum pulse energy required to reach self-injection: 
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where r(l) is the pulse duration after a propagation 
length I. A simple model for the rate of pulse compression 
was put forward in ref. [27] based on the fact that the 
front of the pulse travels at the group velocity of the laser 
in the plasma and the back of the pulse travels in vacuum, 
this produces r(l) « To — (n e l) / (2cn c ). The interaction 
length will be limited by either the length of the plasma 
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FIG. 3. a) electron charge (> 4 MeV) versus aP/P c keep- 
ing the pulse duration constant but varying focal spot qual- 
ity and plasma density (closed circles) or total pulse energy 
and plasma density (open squares) but keeping pulse duration 
constant, b) electron charge versus aP/P c varying pulse du- 
ration while keeping plasma density and energy in focal spot 
constant, c) data from (a) and (b) plotted versus aEn e /n c 
Each data point is an average of five shots and the error bars 
represent one standard deviation. 



target or the pump depletion length l p d ~ cron c /n e 
For the depletion limited case equation [3] reduces to: 
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The threshold density for self-injection for a given experi- 
ment can be calculated from|3]and[4j This model requires 



o 
o 



a 0.5 

CD 








J_ 



J_ 







0.1 0.2 0.3 0.4 

energy within fwhm of spot, aE [J] 

FIG. 4. Observed density threshold as a function of laser 
energy (aE) for our experiment. The solid curve represents 
our threshold model. The dashed curve represents a threshold 
based on aP/P c > 3. 



knowledge of the initial pulse energy, pulse duration and 
the length of the plasma to predict the threshold. As 
equations [3] and [4] are transcendental, the density thresh- 
old for a given laser system must be found numerically. 

A previous study showed that, at low density, the 
threshold is approximately aP/P c > 3 [25], this can be 
rearranged into a similar form to equation |3j 



7re m 2 e c 5 n c 
aE > 3 7> r 



(5) 



We can then use equation [5] to predict the density thresh- 
old for specific experimental conditions. To use this 
model only the initial pulse power is required to calculate 
the threshold density. Combining aP/P c > 3 and equa- 
tion]^ reveals that this threshold model is also equivalent 
to stating that the minimum bubble size for self-trapping 
is constant with density (& p r& > 3.4) in contrast to equa- 
tion CD 

In figure [4] we plot the variation of the observed thresh- 
old density with laser energy (aE). We have defined the 
experimentally observed threshold density as lying in the 
region between the highest density where we observe no 
electron beam and the lowest density where we clearly 
observe a beam. We also show the theoretical threshold 
density based on equations [3] and |4j and the predicted 
threshold based on equation [5] Its agreement with the 
experimental data indicates that our model accurately 
predicts the self-injection threshold, confirming that the 
threshold is reached because the laser pulse undergoes in- 
tensity amplification due to a combination of pulse com- 
pression and self-focusing. 

Our measurements of the threshold density for self- 
injection have been made with only moderate laser pulse 
energies ~ 1 J. Many laser wakefield experiments are 
now being performed with pulse energies ~ 10 J and the 
validity of this model at these higher laser energies can 



5 




S 1 



observed threshold [10 19 cm -3 ] 

FIG. 5. Plot of reported density threshold, n t , versus pre- 
dicted density threshold, n mo dei, for this and other published 
experiments [2l lllll25ll36j . Circles show the predictions of our 
model, Diamonds show the threshold based on aP/P c > 3. 
The line indicates rtt — n mo dei- 



be verified by applying it to previously published data. 
We restrict ourselves to data obtained from experiments 
with gas-jets as guiding structures can affect the trapping 
threshold by changing the way pulse evolution occurs [30] 
or by introducing additional effects such as ionization 
injection [35 j. To calculate the density threshold for a 
particular set of experimental parameters the following 
information is required: the laser energy E, the focal spot 
quality a, the initial pulse duration r and the maximum 
plasma length I. Equations [3] and [4] or equation [5] can 
then be used to calculate the expected density threshold 
for the two models. Kneip et al [2], using a 10 J, 45 fs, 
800 nm laser pulse with a = 0.3, observed a threshold 
density of n e = 2 — 3 x 10 18 cm -3 in an 8.5 mm long 
plasma; our model predicts that the threshold density for 
self-injection should occur at n e « 3 x 10 18 cm -3 . Froula 
et al [25 j, using a 60 fs, 800 nm laser with aE w 6 J 
observed a threshold density of n e w 3 x 10 18 cm -3 in 
an 8.0 mm plasma; our model also predicts n e « 3 x 
10 18 cm -3 . Schmid et al [36] using an 8 fs, 840 nm laser 
with aE ^15 mJ observed electron beams at a density of 
n e ~ 2 x 10 19 cm -3 in a plasma 300 jam long; our model 
predicts a threshold of n e « 2.2 x 10 19 cm -3 . Faure et 
al [11], using a 33 fs, 820 mn laser, reported a dramatic 
decrease in the number of accelerated electrons at n e « 
6 x 10 18 cm -3 in a 3 mm gas jet with aE « 0.5 J. Our 
model predicts a threshold density of n e « 7x 10 18 cm -3 . 

These additional data points, together with those from 
this experiment are presented in figure [5] Due to the fact 
that our model does not depend on a single experimen- 
tal parameter we plot the experimentally observed den- 
sity threshold n t for each experiment on the x-axis and 
against the calculated threshold n mo d e i obtained using 
either equations [3] and [4] or equation [5] Figure [5] shows 
that our model is in good agreement with experiments 



over nearly three orders of magnitude in laser energy, 
whereas the threshold based on equation [5] matches the 
observed threshold over only a very limited range of pulse 
energies: it overestimates the threshold density for low 
energy laser systems and on the other hand would signif- 
icantly underestimate the threshold for very high energy 
laser systems. 

We note that simulations by Yi et al, [22 show that, 
at very low density and an initial laser spot size less 
than the matched spot size, diffraction of the laser pulse 
leads to a lengthening of the bubble which plays a role 
in determining self-injection. In that work they see self- 
injection with a 200 J, 150 fs laser pulse at a density of 
n e = 10 17 cm -3 . Our model predicts that the threshold 
would be n e ~ 4 x 10 1T cm -3 — actually in reasonable 
agreement with [22], however our model relies on pulse 
compression occurring over « 10 cm whereas Yi et al, 
show that in their simulations injection occurs after just 
5 mm . This indicates that our model is only valid for ini- 
tial laser spot sizes greater than or equal to the matched 
spot size (as is the case for the experiments shown in 
figure [5| . 

We now use our model to predict the self-injection 
threshold density for lasers currently under construction. 
For example, our model predicts that a 10 PW laser 
(300 J in 30 fs, A = 0.9 /im, such as the Vulcan 10 PW 
laser at the Rutherford Appleton Lab, or the ELI Beam- 
lines facility in the Czech Republic) could produce elec- 
tron injection at as low as n e ~ 2 x 10 1T cm -3 (assuming 
a = 0.5) in a 6 cm long plasma. For a 1 PW laser (40 J 
in 40 fs, A = 0.8 /im, such as the Berkley Lab Laser Ac- 
celerator, BELLA) our model predicts that self-injection 
will occur at a density of n e ~ 9 x 10 1T cm -3 in 2.4 cm. 

The lower the threshold density of a wakefield accel- 
erator, the higher the maximum beam energy. However, 
for self-injecting accelerators there must be acceleration 
after injection, requiring operation at densities slightly 
above this threshold so that injection occurs earlier in 
the interaction. 

In summary, we have measured the effect of various 
laser parameters on the self-injection threshold in laser 
wakefield accelerators. The simple model we use relies 
on the fact that pulse compression and self-focusing oc- 
cur and that only the energy within the FWHM of the 
focal spot contributes towards driving the plasma wave. 
We find that in cases where the interaction is limited by 
pump depletion, the threshold can be expressed as a ra- 
tio of P/P c , but this ratio is not the same for all laser 
systems: for higher power lasers the threshold occurs at a 
higher value of P/P c than for lower power lasers. When 
the plasma length is shorter than the pump depletion 
length we find that the length of the plasma is an impor- 
tant parameter in determining the injection threshold. 
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